On the analytic structure of scalar glueball operators at Born level 
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We study the analytic structure of the two-point function of the operator F 2 which is expected to 
describe a scalar glueball. The calculation of the involved integrals is complicated by non-analytic 
structures in the integrands, which we take into account properly by identifying cuts generated by 
angular integrals and deforming the contours for the radial integration accordingly. The obtained 
locations of the branch points agree with Cutkosky's cut rules. As input we use different non- 
perturbative Landau gauge gluon propagators with different analytic properties as obtained from 
lattice and functional calculations. All of them violate posit ivity and describe thus gluons absent 
from the asymptotic physical space. The resulting spectral densities for the glueball candidate 
show a cut but no poles for light-like momenta, which can be attributed to the employed Born 
approximation. 

PACS numbers: 11.55.Bq, 02.60.Jh, 12. 38. Aw, 14.70. Dj 



I. INTRODUCTION 

The Standard Model of particle physics describes all 
known elementary particles and their interactions. The 
sector of the strong nuclear force contains quarks and 
gluons, which, however, are not observed directly in ex- 
periments. Searches for partial electric charges, see, for 
example, [H-Q: have all been negative and also gluons 
have not been observed [Bj]. This phenomenon is known 
as confinement and many scenarios trying to explain the 
underlying mechanism exist; for summaries see, for in- 
stance, [6|, LZj. What is observed in nature are bound 
states of quarks and gluons, which are said to be color 
neutral in contrast to their constituents that carry a color 
charge. Composite objects where the constituent degrees 
of freedom are only gluons are called glueballs. They are 
studied by a variety of different approaches, see, e. g., [8] 
for a review. 

Regardless of any details of the confinement mecha- 
nism, one can study the spectral properties of quarks 
and gluons to look for violations of positivity [H, [To| . This 
is considered as a signal of confinement, since it means 
that such a particle has unphysical negative norm contri- 
butions and possesses no Kallen-Lehmann representation 
[ill [l2| . Consequently such particles are absent from the 
physical subspace and are in this sense confined. For 
gluons positivity violations are established from lattice 
[13[ [H| and functional calculations [HI, Qij . It is a non- 
trivial task to obtain physical glueballs from such positiv- 
ity violating constituents. In order to construct meaning- 
ful operators the requirement to have a positive spectral 
density can serve as a guideline. In principle one can ex- 
tract this information by studying the analytic structure 
of the two-point function of such an operator, which was 
done in several studies, see, for example, [l8|. 

The search for bound states with functional approaches 
is usually performed by means of Bethe-Salpeter equa- 
tions, which is a very successful approach for mesons and 



baryons. Only very recently it was extended to gluons 
[l9j . In the present study we take another approach and 
investigate the analytic structure of the two-point func- 
tion of the squared Yang-Mills field strength tensor, the 
F 2 correlator. As a first step we only calculate the ze- 
roth order in the coupling, i. e., we work at the Born 
level. Thus no explicit self-interaction, which would re- 
quire an additional machinery related to renormalization 
[20|, [HI, is considered. However, we account for gluonic 
interactions indirectly by using non- perturb at ive gluon 
propagators as an input which were obtained from the 
interacting theory. Obvious questions are whether the 
signature of a bound state shows up already at the Born 
level in this setting and if there are any major differences 
between the available gluon propagators. 

If the gluon propagator input is simple enough, ana- 
lytic calculations of the spectral representation are possi- 
ble [HSU. However, in general the integrals cannot be 
tackled by hand and a numerical approach is desirable. 
Furthermore, recently numerical solutions for the gluon 
propagator in the complex plane have become available 
[22| which could also be used as an input to numeric cal- 
culations. We thus employ a numerical procedure devel- 
oped specifically for this purpose [23| . The use of Graph- 
ics Processing Units (GPUs) allows a high degree of par- 
allelization and thus a high level of numerical precision. 

Since no generally accepted analytic form of the gluon 
propagator is available, we study several different ver- 
sions. They were obtained from fits to data obtained 
with lattice calculations or Dyson- Schwinger equations. 
Note that they stem from fits of data on the positive 
real axis only and we know from functional analysis that 
this is not enough to fix their analytic structure. First 
numerical results in the complex plane were obtained in 
[22j |. Some fits, however, are motivated by analytic re- 
sults. The most prominent example is the tree- level ex- 
pression of the Refined Gribov-Zwanziger (RGZ) action, 
which can be considered an effective action taking into 
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account the restriction to the first Gribov region and the 
existence of several dimension two condensates |24|, [25| . 
Already the restriction to the first Gribov region leads to 
complex conjugate poles of the tree- level gluon propaga- 
tor [l7j and the condensates modify the gluon dressing 
function further. Since the suggested gluon propagators 
possess quite different analytic structures, it is interesting 
to see how they affect the structure of the F 2 correlator. 
We also study the two-dimensional case. Although not 
physically relevant, it allows us to test the method in a 
simplified setting. 

In Section [TT] we provide the general expression of the 
F 2 correlator in our approximation. We investigate its 
analytic structure in Sections |TTH and [iVl for two and four 
dimensions, respectively, using different gluon propaga- 
tors and conclude in Section IVl 

II. THE CORRELATOR (F 2 (x)F 2 (0)) IN d 
EUCLIDEAN DIMENSIONS 

The quantity of interest is the correlator of the squared 
field strength tensor, which is a candidate for a scalar 



glueball: 

(F 2 (x)F 2 (0)) d = {F^(x)F^(x)F b pr7 (0)F b pa (0)) d . (1) 

d is the (Euclidean) dimension. We will refer to this 
expression as the F 2 correlator. At the Born level 
we neglect the non-Abelian part of the Yang-Mills field 
strength tensor: 

F^ = d»Al - d v A^ (2) 

As customary we will investigate the momentum repre- 
sentation of the correlator: 

(F 2 (x)F 2 (0)) d = J (||^-^( P 2 ), (3) 
where &d(p 2 ) reads [l8| 



/d d k 
^-y d (<?((p - k) 2 )W(k 2 )(k 2 (p - k) 2 + (d- 2){k ■ (p - kj) 2 )) 



(4) 



with C = N 2 — 1 and N c is the number of colors. In the 
following we will use N c = 3 for concreteness but it is 
only a trivial overall factor. Expression (jl]) is valid for 
all kind of gluon propagators which we parametrize as 



p 2 



(5) 



Note that in two dimensions the second term in the in- 
tegrand vanishes. 

In general the integral in eq. (|4]) is UV divergent. To 
render it finite we employ BPHZ renormalization [26l-[28j . 
The highest degree of divergence that appears is of order 
p A . Accordingly we subtract at most up to the fourth 
order in p: 



p=0 



p=0 



(6) 



Note that odd orders drop due to their asymmetric in- 
tegrand. In the following we will drop the superscript r 
again and always refer to renormalized quantities. 

The two-point function A(p 2 ) of a physical operator $ 
(with zero spin) may be written in a spectral representa- 



tion: 



A(p 2 ) 



d d p 
(27r) d 



c ($(;r)$(0)) 



dr 



P(r) 



(7) 



where To is the lowest possible energy of a state. The 
analytic properties of A(p 2 ) are described entirely by this 
expression and can be read off from the spectral density 
p(p 2 ). The spectral density p(p 2 ) must be positive, which 
follows from positivity requirements in quantum theory. 
Violations of positivity are interpreted as an indication 
that such fields are absent from the physical spectrum, 
as found in the_gluon propagator of Landau gauge Yang- 
Mills theory [13I - KL61 122| . The spectral density describes 
the jump across the branch cut along the time-like axis 
and can thus be calculated by 



P(P 2 ) 



1 



lim [A(-p 2 -ie)- A(-p 2 + i e)] (8) 



with — p 2 > tq. Below we will calculate A(p 2 ) for the F 2 
operator, from which we extract the spectral density via 
eq. ©. 
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III. THE ANALYTIC STRUCTURE OF THE F 2 
CORRELATOR IN TWO DIMENSIONS 

We will first test the method in two dimensions. Two 
dimensions are advantageous because renormalization is 
easier: The integrand in eq. (|4]) diverges quadratically in 
p. Furthermore, the second term drops. Two dimensions 
have also the advantage that only one type of solution ex- 
ists [29M3H , whereas in four dimensions two qualitatively 
different kind of solutions can be found, see Sec. IIVI For 
the gluon propagator function &(k 2 ) we use the follow- 
ing fit form motivated originally for four dimensions in 
ref. Ell: 



&(p 2 ) = w\ 
p z 



p A 



• A 2 



1+2k 



(9) 



where k = 0.2 [32[. Note that originally the possibil- 
ity of n = was also considered but can actually be 
excluded because it would require additional unphysical 
prescriptions to be a valid solution |30]. Eq. © describes 
the available data from lattice [33|, [34( and functional 
equations [30j in the IR for w = 1.065 and A = 1 GeV. 
Since deviations from the tree-level for high momenta go 
like 1/(1 + c/p 2 ) the deviations in the UV are negligi- 
ble. In the mid- momentum regime very small deviations 
are observed, which are irrelevant for the present analy- 
sis. However, for easier comparison with d = 4 we used 
here the same parameters as there, viz., w = 2.5 and 
A = 0.51 GeV. 

The results for 0(k 2 ) are shown in figs. [T] and [2j On 
the time-like axis clearly a cut is observed. Since we con- 
sider this as a test case it is important to reproduce the 
expected position for the branch point. From a naive 
guess based on Cutkosky cut rules [35], which, however, 
are not directly applicable here, see, e.g., ref. [36| for de- 
tails, we expect the branch point at p 2 — —1.04 GeV 2 . 
Around this point a structure with four peaks appears, 
which are, however, of finite height as we checked explic- 
itly. A closeup of this region is shown in fig. [3] The 
spectral density, as extracted via eq. (|8j), clearly turns 
negative, see fig. HI There one can also see some numeri- 
cal artifacts, because the discontinuity already starts be- 
low 1.04. However, for larger values, where the numeric 
evaluation is reliable, the spectral density is negative ev- 
erywhere. Thus we conclude that F 2 is not a glueball. Of 
course, this is expected, since in two-dimensional Yang- 
Mills theory no transverse gluons exist which could com- 
pose a glueball and the spectrum must be trivial. 



IV. 



THE ANALYTIC STRUCTURE OF THE F 2 
CORRELATOR IN FOUR DIMENSIONS 



In the four-dimensional case two qualitatively differ- 
ent types of gluon propagators can be found from func- 
tional equations depending on the boundary conditions 
chosen for the equations [16], the so-called decoupling 
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Figure 1: The imaginary part of the F 2 correlator in 2d Lan- 
dau Gauge Yang-Mills theory with the scaling solution gluon 
propagator as input. A branch cut with a singular branch 
point appears on the negative real axis. 
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Figure 2: The real part of the F 2 correlator in 2d Landau 
Gauge Yang-Mills theory with the scaling solution gluon prop- 
agator as input. 



[H, EMI and scaling types solutions 0, 
see also the recent reviews |43|-[45j. The former has a 
non- vanishing value at zero momentum so that £f(0) -1 
can be interpreted as a screening mass, while the latter 
vanishes in the IR as (p 2 ) 2 ^ 1 with k = 0.595353 H2,|46|. 

The tree- level propagator of the Gribov-Zwanziger ac- 
tion reads 



p 



p i + 2N c g 2 ~t 2 ' 



(10) 



where 7 is the so-called Gribov parameter. Qualitatively 
this propagator respects scaling with k — 1. This is true 
also for the ghost propagator, but then also the so-called 
horizon condition must be taken into account which ren- 
ders it then divergent at one- loop level [17[. This was 
confirmed also for two loops [47) . Even more, in a non- 
perturbative treatment the same value for k emerges as 
for the standard Faddeev-Popov action [HI, |49|. Thus 
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Figure 3: A close-up contour plot of the pole structure of the 
singular branch point in the imaginary part of 
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Figure 4: The discontinuity of @(p ) in 2d Landau gauge 
on the time-like axis is negative, thus no Kallen-Lehmann 
representation exists. 



the non-perturbative Gribov-Zwanziger propagators are 
contained in the class of scaling propagators investigated 
here, while the tree- level propagator is qualitatively the 
same but has a different analytic structure. The case 
of the Gribov-Zwanziger propagator can be worked out 
analytically [17j and motivated the introduction of the 
so-called z-particles [l8|, which allow the easy construc- 
tion of an extension of F 2 from the Gribov-Zwanziger 
action that has only physical cuts. The analytic results 
of ref. [l8j were used for testing the presently employed 
numeric code (23| . 

Within the Refined Gribov-Zwanziger framework the 
propagators change their qualitative behavior and be- 
come of the decoupling type [24], |50|. Here it will be 
important that in this case a certain analytic structure is 
determined at tree-level. However, it is not clear if this 
structure is the one of the full non-perturbative propaga- 



tor. The related case of the Gribov-Zwanziger action is 
one example where the difference between tree-level and 
non-perturbative propagator is nicely illustrated [H, |49| . 

For the scaling type we consider the non-perturbative 
part of the gluon propagator suggested in ref. [15] and 
for the decoupling type we investigate the propagator 
of the Refined Gribov-Zwanziger action [24], HBj- The 
latter only provides an effective description in the IR 
and mid- momentum regimes, but the logarithmic run- 
ning of the UV part is not expected to be relevant for 
the present analysis. Furthermore, the non-perturbative 
parts of other suggested fit forms, e. g., in ref. [51[, for the 
decoupling type solution have the same analytic structure 
with only minor changes in the parameters. Thus this 
choice covers a wide range of proposed fit forms. More 
about the specific choices of the fit forms will be said in 
the respective sections below. 



A. Decoupling propagator 

As explained above, one characterizing feature of the 
decoupling solution is a gluon propagator that becomes 
finite at zero momentum. Among possible fits with this 
property a propagator with a constant mass, i. e., 



&(p 2 ) ~ l/(p 2 + m 2 ), 



(11) 



cannot describe the available lattice data satisfacto- 
rily [52], [53j| . A better description is provided by a 
momentum-dependent mass, i. e., m — » m(p 2 ). An of- 
ten employed model for m(p 2 ) is given by [5l| 



m 2 {p 2 ) 



rrin 



p z 



(12) 



where we give only the IR relevant part. Terms taking 
into account the correct UV behavior can be found in 
(5lj , where also the case of an IR constant mass term was 
investigated. Within the RGZ framework the dimension 
two condensates raise the number of possible parameters, 
see refs. [25|, |54( for several possibilities. A fit form mo- 
tivated by the RGZ action that can successfully describe 
lattice data is given by [55| 



&(p 2 ) = C- 



P 



p^ 



u 2 p 2 



■t 2 ' 



(13) 



This propagator describes a more general case than the 
massive propagator in eq. (JTTJ) with the mass model from 
eq. ([T2]) and reduces to that for s = u 2 = t = m 2 ,. For 
the parameters we use the following values from ref. [55| : 
s = 2.508 GeV 2 , t = 0.72 GeV 2 , u = 0.768 GeV and 
C = 0.784. Such a propagator form has also been used 
recently within a toy model including interactions [56| , 
where pole masses could be extracted using a bubble re- 
summation. 

Before we proceed we want to make a few remarks 
concerning the analytic continuation of the momentum 
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Figure 5: The imaginary part of the analytically obtained 
result for the obstructive structure in the complex y-plane for 
x -2 + 2i. 

x = p 2 of the correlator to complex values. As soon as 
x becomes a complex number, the square of the inter- 
nal momentum y = q 2 has to be treated as a complex 
quantity as well. The reason for this is very simple: The 
generic structure of the correlators expressed in hyper- 
spherical coordinates in four Euclidean dimensions is 

J dyy J dz\/l- z 2 f(x,y,z), (14) 

where f{x,y,z) is the regularized integrand whose de- 
nominator can become singular for complex x. z is the 
cosine of the angle between external and internal mo- 
menta. The outer radial integral in y runs along the real 
axis from to some cut-off in the UV denoted by £ 2 . As 
long as x G Mq~, the contour for the radial integral exists 
and the integral can be evaluated. But if x £ Rq~, branch 
cuts induced by the angular integral show up in the com- 
plex y-plane: as the angular integration variable z runs 
through the interval of integration [—1,1], it picks up a 
line of poles for singular values of f(x,y,z). In general, 
the contour along the real axis is obstructed by such a 
branch cut (sometimes more than one) which crosses the 
real axis on the interval [0,£ 2 ]. The contour has to be 
deformed accordingly, in order to avoid these obstructive 
structures. In case of the propagator given in eq. (fT3]h 
the induced obstructive structure in the complex ?/-plane 
is particularly complicated. It features a pair of complex 
conjugate poles as well as two branch cuts. The obstruc- 
tive structure can be worked out analytically by finding 
the poles of the integrand for complex x and y, where z is 
used as a parameter to get the branch cuts. The results 
of the analytic calculation are compared with a numeric 
calculation for the (arbitrarily chosen) point x = —2 + 2i 
in figs. [5]and[6l 

Thus, in addition to deforming the contour around the 
induced branch cuts, the poles have to be treated by ei- 
ther drawing the contours around the poles or by explicit 
subtraction of the residues after the radial integration. 
The latter approach of course also requires an integration 
around the pole in order to obtain the residue. We imple- 




-7 -5 -3 -0.29 1 3 5 7 
Re(y) 



Figure 6: The imaginary part of the radial integrand in the 
complex 2/-plane obtained from numerical integration for x — 
—2 + 2i. Besides the two branch cuts, a pair of complex 
conjugate poles is located at p±^ = —0.29 =b 0.66i. 

mented both approaches and confirmed that they agree. 
However, each method has its strengths and weaknesses 
in certain regions, thus they can be applied according 
to their advantages. Note that the contour depends on 
the specific value of x, since the size and orientation of 
the cuts, but also the residues of the poles depend on 
it. Thus, in order to analytically continue the integrals 
to the complex x-plane, a thorough investigation of the 
complex ?/-plane after the angular integration is neces- 
sary. For more details on the whole procedure and a 
complete worked out example see ref . |23j . 

Further hints of how the resulting structure of the cor- 
relator in the x-plane might look like are provided by the 
Cutkosky cut rules. They allow to calculate the positions 
of the expected branch points. Therefor we need the two 
poles p 2 and p\ of the gluon propagator. They are located 
a t Pi 2 = —0.29 ± 0.66 i. In order to apply the Cutkosky 
cut rules, we take these values to Minkowski space and 
calculate the position of the poles there from 




where i and j are 1 or 2. For i / j we obtain a 
branch point at 2.03. This leads to a cut in Euclidean 
space starting at — 2.03GeV 2 . For i = j the result is 
1.18 d= i 2.63, i.e., there are also two cuts in the complex 
plane starting at —1.18 =F z2.63GeV 2 . These cuts are 
considered unphysical in the sense that they forbid to 
write a spectral representation for the correlation func- 
tion of F 2 . One motivation for the introduction of i- 
particles in ref. [18| was that they provided an easy way 
to get rid of such unphysical cuts. If these cuts were 
absent, we would have a proper spectral representation 
for the glueball, since the discontinuity across the cut on 
the time-like axis leads to a positive spectral density as 
shown in fig. [T2J The locations of the branch points as 
predicted from the Cutkosky analysis is also manifest in 
the obstructive structure in the complex ?/-plane. Upon 



6 



6 



4 




-4 



-6 

-6 -4 -2 2 4 6 
Re(y) 

Figure 7: At x = — 2, the endpoints of the cuts in the y- 
plane are "touching" the poles. A valid contour connecting 
to the UV cut-off £ 2 must in one way or another be drawn 
through one of these points. This restriction leads to a non- 
analyticity. In fact, this point has also been obtained from a 
Cutkosky analysis. 
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Figure 8: At x — —1.18 — 2.72 a valid contour must be drawn 
through the pinching curly cut, whose one end-point coincides 
with the pole. Again, the empirically determined value of x 
where this happens is close to the Cutkosky prediction for the 
branch point. 

varying the complex values of x we could find three values 
of x for which the obstructive structure in the ?/-plane is 
such that it leads to a non-analyticity. They are depicted 
in figs. [7J [8] and [9] The values were determined empiri- 
cally and are in agreement with the three values obtained 
from the Cutkosky analysis. The empirical procedure of 
finding the pinching points gives rise to a small deviation 
from Cutkosky's predicted values. If they could be de- 
termined with arbitrary precision, they would coincide. 
Non-trivial analytical properties at these points are ex- 
pected because in all three cases the integration contour 
must be in one way or another drawn such, that it runs 
through the endpoint of one of the cuts which happens to 
lie just on top of one of the singularities. A possible path 
connecting to the UV-cutoff £ 2 has to pass this point, 
which leads to a non-analyticity. Thus, the three points 
can also be found by looking for values of x for which the 
possible contours have to be drawn exactly through the 
point where one (non-singular) branch point "touches" 
one of the poles. 

In figs. [lOj [11] and [12] the numerical results for the 
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Figure 9: For x — — 1.18 + 2. 7i the situation is similar to the 
one found for the complex conjugate of x, as depicted in fig. [8] 



imaginary part, the real part and the discontinuity along 
the physical cut of the correlator are shown, respectively. 
Because of the very restrictive analytic structure in the 
complex plane of the radial integration variable the ex- 
act structure is numerically almost not resolvable. Both 
methods described above give the same results in well- 
behaved regions, but tend to get unstable exactly in the 
directions of the three cuts. For the unphysical cuts 
this is simply due to the fact that the positions of the 
poles induced by the propagators have the same argu- 
ment as the branch points. As the branch cuts always 
open into the direction given by arg(x), the only way the 
contour can be drawn out of the cuts is by going into 
the direction given by arg(x), i.e. running straight into 
the poles. All points having the same argument as the 
pole-positions suffer from this problem. The unphysical 
cuts also have the property that they open very slowly, 
which makes a precise numerical determination of the 
branch point locations almost impossible. Nevertheless 
the small, but non-vanishing discontinuity clearly indi- 
cates the presence of the unphysical cuts. For small nega- 
tive real numbers and (almost) vanishing imaginary part, 
the contour runs very close to the obstructive branch cut 
over large distances, which makes it also very hard to 
obtain precise values in this region. Because of this, the 
physical branch cut opens too early as compared to the 
Cutkosky prediction. This stability issues can also be 
seen in the discontinuity of the physical cut, as shown in 
fig. [12] The few points below the axis can be attributed 
to the unstable numerics, thus the spectral function is 
strictly positive. The physical cut is clearly visible in 
fig. [10] We want to stress that the results can be im- 
proved by adjusting the integration paths in the ?/-plane. 
However, this complicates the integration routine further 
and requires considerable fine-tuning and more integra- 
tion path classes. To illustrate this we compare the effort 
required for the i-particles calculation of [23| with that 
of the Refined Gribov-Zwanziger case: For the former, 
three different types of integration paths are necessary in 
order to avoid the singularities in the ?/-plane. Each of 
them consists of two to five integrations. For the Refined 
Gribov-Zwanziger propagator we employed here five dif- 
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Figure 10: The imaginary part of the F 2 correlator in 4d 
with Refined Gribov Zwanziger gluon propagators. The un- 
physical cuts open very slowly, making it almost impossible 
to determine the branch point location numerically. 
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Figure 11: The real part of the F 2 correlator in 4d with Re- 
fined Gribov Zwanziger gluon propagators. 

ferent classes of paths. In order to obtain better results 
we expect that at least two more are necessary. Since it is 
not evident that the fit given in eq. ([T3]) is the true non- 
perturbative form, as indicated also by numeric results 
for complex momenta [22], we refrain from this further 
complication. Nevertheless, we could demonstrate how 
the cuts in this case arise and obtained results for the 
branch point positions in agreement with the Cutkosky 
cut rules. 



B. Scaling propagator 

The non-perturbative part of the four-dimensional 
scaling gluon propagator can be described by [HI 




1 2 3 4 5 
-p 2 [GeV 2 ] 

Figure 12: The discontinuity of 0(p 2 ) for the 4d Refined Gri- 
bov Zwanziger propagators. Some points are below the axis 
because of stability issues of the numerics. For the same rea- 
son the cut opens way too early as compared to the Cutkosky 
prediction. 

where n — 0.595353 [H, |46|. In order to account for 
the full momentum dependence of the gluon propagator, 
this function has to be multiplied by another function 
that describes the UV logarithmic tail [15]. However, we 
are here only interested in the non-perturbative content 
of the propagator and thus use expression ([T6]) . This 
also avoids technical complications related to the defor- 
mation of the contour in the integration process. For 
the parameters we use the value from the fit in [15], viz. 
A = 0.51 GeV and w = 2.5. 

The resulting correlation function of F 2 is shown in 
figs. [T3l and fl4l As expected, there is a branch cut on the 
time-like axis, see fig. [13l The jump of the correlation 
function, which is directly related to its spectral density, 
is plotted in fig.[15j In this case the numerics are perfectly 
under control and we can easily deduce the position of the 
branch point, which is 1.04 GeV 2 . We thus confirm that 
the F 2 operator corresponds to a physical quantity with 
a positive spectral density. However, we do not observe 
any poles which, of course, are not expected to appear at 
Born level. 



V. SUMMARY AND CONCLUSIONS 

In this work we investigated the analytic properties of 
the correlation function of the operator F 2 , which is a 
candidate for a scalar glueball. In our setup we used 
different non-perturbative gluon propagators which are 
obtained from fits to lattice and/or continuum results. 
For some simple cases the calculations can be done an- 
alytically. However, for future investigations numerical 
methods are required for several reasons. First of all, 
one would like to directly use results for the gluonprop- 
agator. Such data became available only recently [22|. In 
this case no fits are available and numeric data is the in- 
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Figure 13: The imaginary part of the F 2 correlator in 4d 
Landau Gauge Yang-Mills theory with scaling solution gluons 
as input. 
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Figure 14: The real part of the F 2 correlator in 4d Landau 
Gauge Yang-Mills theory with scaling solution gluons as in- 
put. 



put. Second, extensions of the current scheme most likely 
can only be calculated numerically. With the present 
work we have laid the basis for such calculations. 

A new complication, arising in future calculations 
when interactions are taken into account, is the neces- 
sity of renormalization. The required machinery for this 
is mostly available [20|, [2l| , but calculations will certainly 
become more complicated. In this respect it will be use- 
ful to continue as done here with GPUs as they allow 
a great deal of parallelization, which in this case can be 
done straightforwardly as no communication between dif- 
ferent threads is necessary. 

Since we worked at the Born level, we did not find 



any poles but only cuts for the F 2 correlator. In case of 
the scaling type gluon propagator there was only one cut 
on the time-like axis. For the decoupling type propagator 
two additional cuts in the complex plane were found that 
forbid a spectral representation. However, at this stage 
we can not draw any physical conclusions from this differ- 
ence between the two solutions, because the unphysical 
cuts can directly be attributed to the analytic properties 
100 




Figure 15: The discontinuity of &(p 2 ) using scaling solution 
gluon propagators. 



of the decoupling fit. While the employed fit form can de- 
scribe lattice data reasonably well, this seems no longer 
to be true for complex momenta. Indications for this are 
also found by a direct calculation of the gluon propaga- 
tor (22[, where no poles in the complex plane have been 
found. 

Several possibilities for future calculations exist. On 
the technical side we have the inclusion of interactions, 
while on the formal level a better input is required for 
the gluon propagators. Due to recent results [22[ this 
may even come from calculations directly for complex 
momenta instead of fits to Euclidean momenta. The work 
presented here will in any case provide useful guidance 
for determining the analytic structure of correlators. 
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